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ABSTRACT : Isoscalar dipole (compression) mode is studied first using schematic harmonic- 
oscillator model and, then, the self-consistent Hartree-Fock (HP) and random phase ap- 
proximation (RPA) solved in coordinate space. Taking 208p]^ g^j^^ ^j^g SkM* interaction 
as a numerical example, the spurious component and the strength in the continuum are 
carefully examined using the sum rules. It is pointed out that in the continuum calcula- 
tion one has to use an extremely fine radial mesh in HP and RPA in order to separate, 
with good accuracy, the spurious component from intrinsic excitations. 



PACS numbers : 21.10.Re, 21.60.Jz, 23.20.Js, 27.80.+W 



1 Introduction 



Though the isoscalar (IS) dipole resonance (a compression mode) was theoretically studied 
already more than 20 years ago H, ^, ^ and various experimental efforts to pin down the 
IS giant dipole resonance (ISGDR) have been made [|, ||, ||, |I0|, the distribution of the 
IS dipole strength obtained from the analysis of experimental data is still under dispute. 
The difficulty in extracting the experimental dipole strength comes from the ambiguity 
in the parameters of the optical potentials ||TT|, which are needed in the analysis of the 
relevant {a, a) scattering data. In order to simplify our discussion, as a numerical example 
in the present article we discuss only the doubly magic nucleus g2^Pbi26. A current hot 
issue in both theory and experiments is whether or not a considerable amount of IS dipole 
strength is found in the energy region of £^=8—17 MeV, which is much lower than the 
so-called ISGDR with a peak around £'^=23 MeV. Available calculations have a tendency 
to give 1-3 MeV higher peak-energy of ISGDR than that identified so far in experiments. 

The calculation of IS dipole strength requires special care since the spurious translation 
mode has to be fully separated from intrinsic excitations. In the self-consistent random- 
phase-approximation (RPA) calculation the spurious mode should appear at zero energy 
with a finite contribution to the energy- weighted sum-rule (EWSR). The response function 
to the IS dipole operator 

D = j:irlYiM) (1) 

k 

contains the strengths of both the spurious mode and intrinsic excitations. In actual 
numerical calculations of self-consistent RPA the spurious mode appears not at zero energy 
but at either a finite energy or an imaginary energy, due to numerical inaccuracy as well 
as some possible practical approximations. This may lead to a question of whether or 
not the calculated IS dipole strength obtained in the energy region of intrinsic excitations 
is free from the spurious component. The difficulty in the present problem lies in the 
fact that a very small amount of the spurious component remaining in the energy region 
of intrinsic excitations can in turn produce an appreciable amount of the strength for 



the operator (0) due to the power of r^. If the transition density coming from spurious 
component remaining in the region is proportional to ^ where po expresses the ground- 
state density, the IS dipole strength which is free from the spurious component is obtained 
by using the operator 

D = J2irl-vr,)Y,o{9k) (2) 

k 

where ri=^{r'^) ^ instead of using that in (|l]). However, if some spurious com- 
ponents remain at the energies of various unperturbed particle-hole (ph) excitations, the 
relevant transition densities are similar to those of respective original ph excitations rather 
than If so, using the operator in (H) may not help to take away all strength coming 
from the spurious component remaining in the region of intrinsic excitations. 



In ref. [jT2[ we studied both the IS and isovector (IV) dipole mode using the self- 
consistent RPA calculation with the SkM* interaction, in which both the IS and IV 
correlations in RPA are taken into account solving both the Hartree-Fock (HF) and 
RPA equations in coordinate space using the Green's function. Examining the IS dipole 
strength obtained by using the operator in (|^), in the energy region of 7 < Ex < 17 MeV 
we obtained an appreciable amount of IS dipole strength almost exactly at the energies 
of various unperturbed ph excitations. Since the residual interaction relevant to the IS 
dipole correlation is in general attractive, the strength increased by the RPA correlation 
is expected to be found at energies lower than the unperturbed ones. Therefore, using the 



procedure described in ref. [I^ we tried to remove the strength at finite energies, which 
seemed to come from the remaining spurious component. In the continuum at a given 
energy the excitations of various modes with various quantum numbers may coexist. The 
separation of the strengths coming from modes with different quantum numbers is easy, 
however, we have found no established way of separating the strengths originating from 
different modes with the same quantum numbers. Since both the spurious mode and IS 
dipole resonance are isoscalar and have the same spin-parity, f = 1~, our basic idea used 



in ref. was the following. Near the isolated resonance, E ^ En, the relation 

Im {GnPAir, P, E)) oc p'^ir) p'^i?) (3) 

is known. At a high energy in the continuum, where resonances i and j with respective 
finite widths may overlap, we may write 

Im[GRPAif,r',E)) oc a, pf{f) pf {?) + b, 6pfir) 6pfir') 

= F{f,r') (4) 

where the modes, i and j, have the same spin-parity and are diagonal in RPA. Choosing 
some /(r) so that 

' dr' pf{r') f{r') = (5) 



one obtains 

5pf (f) oc / dr' F{r,r') f{r') (6) 



Choosing /(r) = rYi^ in the IS dipole case, we obtained Spf^f), which was supposed to 
be the transition density at the energy E, that carried the spurious rYi^ strength. Then, 
we subtracted Spf{r) from the transition density ptr{r) at E for the operator D in ([|), 

pt,(f) oc J d? Im (GRPA{f, f', E)) {r'f (7) 



Using the above subtraction procedure, it seemed that in ref. |]T2[ we subtracted a bit too 
much strength as was criticized, for example, in ref. ||T3[. The oversubtraction happened, 
since the remaining ph excitations which carry some spurious strength were not really 
RPA solutions. 

In all available calculations of ISGDR except those in refs. |12|, |l^ the states in 



the continuum are approximated using the expansion in terms of discrete basis. Namely, 
the continuum wave functions are expanded in terms of harmonic oscillator basis with 
a maximum principal-quantum- number N^ax- In those discrete-basis calculations the 
operator in (0) is commonly used in order to subtract the spurious components. The IS 



dipole strength obtained in the energy region of i^^ = 8 — 17 MeV depends very much 
on pubhcations. The relativistic RPA calculation in ref. [0 produced especially a large 
amount of the IS dipole strength in the low-energy region, for all effective interactions 
used. In none of the discrete-basis calculations the EWSR carried by the spurious state 
and intrinsic excitations is carefully examined. 

In the present paper we use EWSR as a measure of the validity of numerical calcula- 
tions. The EWSR for the IS dipole operator in (|l]) is written as 

n 

where En expresses the energy of RPA states, while {r'^) on the r.h.s. is evaluated with 



the HP ground state |14|. The total EWSR in (§) consists of the contribution by the 
spurious state Sgpr and that by the intrinsic excitations Sintr, S = Sspr + Sintr- The latter 
can be written as [^] 

intrinsic,i 

3 ^ f fl A \ ( ^ ^ , ^> 2\2 



, > - — ) (9) 

which is indeed equal to the EWSR for the IS dipole operator in (0). 

In sect. 2 we use a harmonic oscillator model to evaluate the distribution of the 
EWSR of IS dipole strength among the spurious state, the lower-lying and higher-lying 
intrinsic excitations. This schematic model is very much simplified, nevertheless, we 
expect that the semi-quantitative feature of the result obtained from the model should 
remain unchanged in realistic calculations. In sect. 3 we present the result of the self- 
consistent RPA calculation with the Skyrme SkM* interaction on ^°'^Pb, which is the 
continuum calculation solving both HP and RPA in coordinate space. The result is 
compared with discrete-basis calculations, in which the wave functions in the continuum 
are expanded in terms of harmonic oscillator basis. In sect. 4 conclusion and discussions 
are given. 



2 Harmonic oscillator model 



The distribution of IS dipole strength may be studied in terms of a simphfied model |15 



in which the particle motion is described by a harmonic oscillator potential. The particle 
excitations produced by the field (0) are governed by the selection rule AA^=1 or 3 and 
have the energies hujQ and 3huJo, respectively. For a single completed shell with principal 
quantum number N, the transition strength is given by 

E \{^N+i\r'Y,oWN)\' = ^ + 1){N + 2){N + 3)(13iV2 + 52N + 50) 

(10) 



E \{^N+3\r''YwWN)\^ = ^ I: 



h V 1 



— {N + l){N + 2){N + 3)iN + A){N + 5) (11) 



2MuJoJ 12 

where represents the quantum numbers needed to specify the single-particle states in 
the shell A^. For simplicity, we do not include the spin-isospin degeneracy in the present 
model. Since the last three occupied shells contribute to the total AA^=3 strength we 
obtain 

Nf Q / ft \ ^ 

E E 

N=Nf-2 l^N,VN+i 



Nf 3 f n V 1 

E E IK+3k'>^iok^)P = -hj7l— -{Nf+1){Np+2){Nf+3){3N'p+12Nf+20) 
=Nf-2 un,um+3 47r \ZMuJoJ LZ 



(12) 

while only the last-filled = Np shell contributes to the AA^=1 strength. 

Now, the major part of the AA^=1 strength, which is obtained by setting = Np in 
(p!OD , corresponds to the spurious excitation, namely center of mass motion, while in the 
present harmonic oscillator model the AA^=3 strength contains no spurious component. 
Subtracting the strength coming from the spurious excitation 



n V25 



— {Np + l){Np + 2Y{Np + 3), (13) 



471 \2Muo) 24 

from the expression (ITUp with = Np^ we obtain the strength of the intrinsic AA^=1 
excitation 

^ / f) \ ^ 1 

-Np{Nf + 1){Nf + 2){Np + 3){Np + 4) (14) 



47r \2Muoq 24 



We note that the total number of particles in the present model is written as 



1 



{Nf + 1){Nf + 2){Nf + 3) 



(15) 



6 



while 



N=0 



( 



2MuJol 4 



iNF + l){NF + 2f{NF + 3) 



(16) 



and 



E = TTTT- + + 2)(^i^ + 3)(iV| + 4iV^ + 5) (17) 




The spurious contribution ([T3|) is estimated by recognizing that the spurious compo- 



nent is found only at the AA^=1 excitations and that the transition density of the spurious 
excitation is proportional to the radial derivative of the ground-state density, dpo/dr. 

From the expressions, ([I^ ) and ([T^) , one obtains the ratio of the unperturbed transition 
strength of AA^=1 excitations to that of AA^=3 ones to be 1:6 for Np ^ 1. Correspond- 
ingly, the ratio of the unperturbed EWSR is 1:18. 

Introducing the residual interaction of a separable form 



where c expresses an attractive coupling constant, we solve the RPA equation and estimate 
the distribution of the transition strength together with the RPA energies. The following 
results of cases (a) and (b) are summarized in Table I. 

(a) If we choose c so that the higher-lying RPA solution is obtained at {2.5)hujQ, the 
lower-lying RPA solution is found at (0.872)^Ci;o. Then, the RPA transition strength is 
distributed with 1.09:2.79 over the lower-lying and higher-lying RPA solutions, while the 
distribution of the EWSR is 0.95:6.98 . 

(b) If we choose c so that the higher-lying RPA solution is obtained at (2.8)^a;o, the 
lower-lying RPA solution is found at (0.961)/;,ti;o. Then, the RPA transition strength is 
distributed with 0.60:2.62 between the lower-lying and higher-lying RPA solutions, while 
the distribution of the EWSR is 0.58:7.34 . 



V = cD^D 



(18) 



Comparing the peak energies of the RPA solutions with those obtained by available 
self-consistent HF plus RPA calculations, one expects that the realistic situation may lie 
somewhere between the case (a) and (b). Before the spurious component 6 pf{r) in (|^) was 



subtracted in our self-consistent RPA calculations of ref. |T2|, a number of lower-energy 
peaks at Ex < 17 MeV for the operator were found nearly at respective unperturbed ph 
energies. In contrast, the higher-energy peak around 25 MeV was collectively constructed 
by shifting the unperturbed strength coherently from the region of higher energy. In 



other available self-consistent RPA calculations, for example in refs. [0 and [17 



m 



which the expansion in terms of harmonic oscillator basis was used, or in ref. [|I^ which 
was a continuum calculation of relativistic Hartree plus RPA calculation, a considerable 
amount of low-energy strength appeared also around unperturbed ph energies. Compared 
with those HF plus RPA calculations, in the case (a) of Table I the higher-lying RPA peak 
is pushed down by about the same order of magnitude and the lower-lying peak is too 
strongly pushed down, while in (b) the higher-lying peak is appreciably less pushed down 
and the lower-lying peak is still slightly lower. 

On the other hand, one might wonder whether or not the ratio of the coupling between 
the AA^=1 and 3 excitations to that within the excitations of a given AA^ is realistic in 
the present schematic interaction of (0). In the present harmonic-oscillator model the 
effective strength of the coupling between the AA^=1 and 3 excitations can be measured 
from the increase of EWSR in the lower-lying RPA solution compared with that in the 
unperturbed AA^=1 excitation. The portion of the EWSR in carried by the lower- 
energy peak is increased from 2.3 to 5.2 percent in the case (a) and to 3.2 percent in the 
case (b), as shown in Table I. 



3 Hartree-Fock and RPA calculations with Skyrme 
interactions 



In the HF calculation with the SkM* interaction for ^osp]^ obtain 

(r^) = 30.85 fm^ 

(r^) = 1218.2 /m^ (19) 

Thus, using eqs. (Q) and (H) the EWSR becomes 

S{r^Yio) = 13.79 x 10^ fm^ MeV (20) 

Sintrir^YiQ) = 5.632 x 10*^ fm^ MeV (21) 
From eqs. (pOj) and (^Tf) the contribution by the spurious state to the EWSR for the 



operator r^Yi^ is 8.16 x 10^ fm^ MeV, which is 59.2 percent of the total EWSR in (|0D. 

In the present work we perform the self-consistent RPA calculation with the SkM* 
interaction, in which both the IS and IV correlations in RPA are taken into account 
solving both HF and RPA equations in coordinate space using the Green's function. The 
difference of the present calculation from that in ref. []T2[ is that we use the radial mesh 



(Ar) of 0.1 fm in both the HF and RPA calculations and, moreover, we carefully check the 
sum-rules consumed by both intrinsic and spurious excitations. We present the numerical 



result without adopting the subtraction procedure, which was used in ref. [|T2[ and was 
also explained in the Introduction. In ref. |]12| Ar=0.1 fm was used in HF, while we 
adopted Ar=0.3 fm in RPA. In Fig. 1 we show the RPA response functions, which are 
smeared out using the width of 0.5 MeV, for the operator (|I|) obtained with Ar=0.1 fm 
and 0.3 fm and for the operator (^ calculated with Ar=0.1 fm. A great, somewhat 
unexpected improvement in the result of the present calculation with a finer radial mesh 
in RPA is : the appreciable difference between the response functions to the operators (|l]) 
and @ in the energy region of the higher peak around 25 MeV is drastically decreased. 
Moreover, the spurious state appears now slightly below MeV (namely, at an imaginary 
energy) and consumes a much larger portion of EWSR, Sspr- 



As in ref. |]T2l, our effective interaction used in RPA consists only of the part in tfie 
spin-independent channels. Using discrete-basis calculations, in which the inclusion of 
the full Skyrme effective interaction is straightforward and easy, we have checked that the 
RPA response function of the present IS dipole operator changes very little when we use 
the full effective interaction, as seen in Tables II and III. 

In Table II we show the transition strength for the operator (|l]), the contribution to 
EWSR, and the percent of the total EWSR, S'(r^Y'io), at the spurious state and in the 
region of 2 < E^. < 17 MeV, 17 < < 70 MeV, and > 70 MeV. In Table III the 
same quantities are tabulated for the operator (|2|). The borderline 17 MeV is chosen, 
since almost all AA^ = 1,1'^ = 1^ ph excitations in the HF potential lie below 17 MeV 
and no unperturbed 1~ ph excitations are found for 18 < -E^ < 24 MeV. 

As is previously known, it is seen from Table II that the portion of the EWSR found 
in the energy region of intrinsic excitations {Er^ > 2 MeV) clearly shows the presence of 
the spurious components remaining in the region. From Table III it is seen that the ratio 
of the EWSR lying in the region of E^ < 17 MeV to that for E^ > 17 MeV is 1:4.9. More 
importantly, the portion of the EWSR obtained for E^ > 2 MeV is 39.1 percent of the 
total EWSR, compared with the full value, 40.8 percent. It should be mentioned that in 
the continuum calculation it is technically very difficult to obtain accurate values of the 
transition strength carried by very sharp peaks including the spurious state, in contrast 
to the discrete-basis calculations. Furthermore, in the present continuum calculation the 
spurious state appears slightly below MeV. Thus, reducing the effective ph interactions 
by multiplying a factor of 0.9902, we obtain the spurious state at i?2:=0.107 MeV, from 
which we estimate the carried transition strength. Though it is clear that the continuum 
calculations performed using a much finer radial mesh than 0.1 fm would further improve 
the accuracy, here we do not pursue it any further. 

In Fig. 2 we show the calculated RPA response functions to the operators, (|l|) and 
(H), as a function of excitation energy. Though the peaks at E^ > 7.5 MeV lie in the 



continuum, calculated lower-lying peaks are too sharp to be plotted. Since the spreading 
width is anyhow not included in the RPA response function, in Fig. 2 we plot the 
quantities, which are obtained by smearing out the calculated RPA response function 
with the width of 1 MeV. It is seen that an appreciable amount of IS dipole strength is 
present in the energy region of Ex=7-17 MeV, which is much lower than the so-called 
ISGDR around 25 MeV. The presence of the lower-lying IS dipole strength indicates that 
the ISGDR appearing around 25 MeV is not really a good collective mode in the nuclear 
system. 

For reference, in Tables II and III we show also the result of discrete-basis calculations 
for the operators, (|Tp and (|^). How to truncate the space of harmonic oscillator basis seems 
to depend on publications. What we have adopted in the present work, which is along the 
line of ref. [|^, is the following. First we construct the HF potential, which is calculated in 
coordinate space and, then, calculate the one-particle energies and wave-functions for the 
potential using the harmonic oscillator basis with the radial node n^'^ < n^^^, where n^^^ 
is taken to be independent of £. Then, the RPA equation is solved including the lowest- 
lying \{N^^^ —i) one-particle levels for a given i, where A^^ox is taken to be independent 
of i. This truncation of the ph space included in RPA is approximately the same as that 
limiting the ph space by the ^-independent maximum energy of particle configurations. 
It is noted that, as n^^^. increases, the bound one-particle energies and wave-functions 
approach the eigen energies and eigen functions of the HF potential, respectively, while 
the properties of the positive-energy one-particle levels can be far away from the correct 
ones unless ^i^oL is taken to be cxd. In actual calculations one must take finite values of 
'^moL ^max ^^^1; thus, it is uot guaranteed that using finite larger values of and 
^max produces a more accurate result of IS dipole strength. 

For the parameter set, n^^c = 12 and A'mox = 24, in Tables II and HI we show both 
the result obtained by using the full effective interaction in RPA and the one by employing 
the same effective interaction used in our continuum RPA calculation. The fact that very 



little difference is found between these two results indicates that the IS dipole strength 
obtained from our continuum RPA calculation is reliable. Comparing the numbers in 
the last column of Table II with those in Table III, in the discrete-basis calculation the 
difference between the strengths calculated with the operators (|l|) and (^) is relatively 
minor in the region of > 17 MeV. On the other hand, from Table II we notice that 
the calculated energy of the spurious state is pretty high compared with the continuum 
calculation and, furthermore, only 92 percent of the correct portion of EWSR is consumed 
by the calculated spurious state. The properties of the spurious state converge surprisingly 
slowly to the right ones, as we increase n^^c N^^^. On the other hand, if we use a 
value of larger than 12, the calculated response function in the continuum starts to 
behave strangely instead of converging to some reasonable result. This strange behavior 



of discrete-basis calculations seems to be known for experts in the fields ||19 |. 

4 Conclusion and discussions 

First, the distribution of IS dipole strength for heavy nuclei is studied using the schematic 
harmonic-oscillator model, in which the strength is distributed among the spurious mode, 
the low-energy and high-energy peaks. Though 96 percent of the unperturbed AA^ = 1 
strength for the operator r'^Yi^ corresponds to the spurious excitation, the remaining 4 
percent contributes to the intrinsic excitations. The spurious state consumes 56.8 percent 
of the total EWSR for the operator r'^Yi^, while the low-energy RPA peak may carry up 
till 5 percent of the EWSR. 

Secondly, the result of IS dipole strength in ^°^Pb is presented, which is calculated us- 
ing the self-consistent HP plus RPA solved in coordinate space using the Green's function 
method, with better accuracy than in ref. [|l^]. If one performs numerical calculations 
with very good accuracy, the IS dipole strength of intrinsic excitations should not depend 
on whether the operator (|l|) or is used. At the same time, the spurious mode should 
appear close to zero energy and consume the whole portion of EWSR expected. In our 



present continuum calculation we could not attain such a satisfactory accuracy. On the 
other hand, if the calculated strengths with the operators and @ are appreciably 
different from each other in the energy region of intrinsic excitations, the spurious compo- 
nent remaining in the region may not have the transition density proportional to In 
the present work we have not attempted to take away this kind of possible spurious com- 
ponent from the region of intrinsic excitations, since we find no reliable way of excluding 
it. 

Using the radial mesh of 0.1 fm in both the HF and RPA calculations, the spurious 
state is obtained slightly below MeV. Multiplying the effective ph interaction by a factor 
of 0.9902, we push up the energy to 0.107 MeV and estimate the EWSR carried by the 
state, which turns out to be the major portion of the expected EWSR. On the other 
hand, the EWSR carried by the calculated states for > 2 MeV is 114 and 96 percent 
of the expected EWSR for the operators (|l|) and (H), respectively. The accuracy of the 
continuum calculations is appreciably improved by decreasing Ar used in solving RPA 



from 0.3 fm in ref. |12| to 0.1 fm, however, it is still away from being very satisfactory. 
In the continuum calculations the accuracy of separating the spurious component from 
intrinsic excitations is surprisingly sensitive to the size of radial mesh employed in RPA. 
Thus, the response function of IS dipole strength has to be calculated using a very fine 
radial mesh, in order to obtain a reliable numerical result. From the present study we find 
that slightly more than 5 percent of the total EWSR is expected as intrinsic excitations 
lying in the low-energy region below < 17 MeV. 

The expansion of HF and RPA solutions in terms of harmonic oscillator basis can 



work for bound states if one includes a sufficiently large number of bases [^. Thus, 
in the present study of IS dipole strength the properties of the spurious mode become 
better and better as the number of bases increases, though we find that the convergence 
of the properties is unexpectedly slow as n^^c ^max increase. On the portion of 
EWSR found in the continuum the discrete-basis calculation presented in Table II and 



Fig. 1 produces 9 percent larger in the 2 < < 17 MeV region and 6 percent smaller 
for Ex > 17 MeV, compared with our continuum calculation. However, the real problem 
is : If we further enlarge the space of discrete-basis, the calculation starts to produce 
more than 100 percent of EWSR, which are certainly unreasonable. Thus, it is not clear 
whether or not the discrete-basis calculations can be reliably used for the estimate of IS 
dipole strength. 
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Figure captions 

Figure 1 : RPA response functions to the IS dipole operators, (|l]) and (H), as a function of 
excitation energy, which are obtained from the self-consistent HF plus RPA calcula- 
tions solved in coordinate space. A radial mesh Ar=0.1 fm is used in HF, while the 
results obtained by using Ar=0.1 and 0.3 fm in RPA are compared. We show the 
calculated response functions, which are smeared out using the width of 0.5 MeV. 

Figure 2 : RPA response functions to the IS dipole operators, (|1|) and (0), as a function of 
excitation energy, which are obtained from the self-consistent continuum calculation 
(CRPA) with Ar=0.1 fm and the discrete-basis calculation with n^^^ = 12 and 
N^^^ = 24 (DRPA). In all calculations the (1, r ■ t) effective interaction (see the 
caption to Table II) is used. In order to see more easily the total distribution of the 
strength, we plot both CRPA and DRPA response functions, which are smeared out 
using the width of 1 MeV. 



Table 1: Energy, reduced transition strength B{D), and contribution to the EWSR for 
the operator r^Yi^, which are calculated in the harmonic oscillator model for Np ^ 1. 
B(D) and EWSR are given in arbitrary unit. 



Case 


energy in huJo 


B(D) 


EWSR (%) 


unperturbed 


spurious (l.U) 
1.0 
3.0 


^ = 1.042 

T4 = 0-042 
1 = 0.25 


^ = 1.042 (56.8) 
2^ = 0.042 (2.3) 
f = 0.75 (40.9) 


(a) 


0.0 
0.872 
2.5 


oo 
0.109 
0.279 


1.042 (56.8) 
0.095 (5.2) 
0.698 (38.0) 


(b) 


0.0 
0.961 
2.8 


oo 
0.060 
0.262 


1.042 (56.8) 
0.058 (3.2) 
0.734 (40.0) 



Table 2: Energy, reduced transition strength B{D), and contribution to the EWSR for 
the operator r^Yi^ in ^°^Pb, which are calculated in the self-consistent RPA using the 
SkM* interaction. Inside the bracket of the forth column the ratio to the total EWSR for 
the operator r^Fi^, (pO]) , is expressed in percent. For discrete calculations with n^^^ = 12 
and N^^^ = 24 we show two sets of calculated results, the one obtained by using the full 
effective interaction in RPA and the other calculated by using only the spin-independent 
part as that employed in our continuum RPA calculation. The former is denoted by 
"full", while the latter is marked by "(1, r ■ r)" . In the continuum calculation it is 
technically very difficult to estimate an accurate value of the transition strength carried 
by the spurious state, which is given with in the third and fourth column. See the 
text for details. 



Type of calculation 


energy in MeV 


B(D) in fm^ 


EWSR in fm^MeV (%) 


correct calculation 




0.0 


oo 


8.16X 10^ (59.2) 








< 




5.63X 10*^ (40.8) 




continuum 




0.107 


^ 7.5x 10^ 


^ 8.0x 10^ 




(l,r-r) 




2<E^< 17 


1.275 X 10^ 


1.433x10^ (10.4) 








n <E^< 70 


1.973 X 10^ 


4.856X 10*^ (35.2) 
1.14 X 10^ (0.83) 








70 < < 150 


1.3 X 10^ 


discrete 


\ ^max y max 


= 16) 


1.401 


5.342 X 10^ 


7.485 X 10^ (54.7) 




full 




2 < < 17 


1.347 X 10^ 


1.534 X 10^ (11.1) 








n <E^< 70 


1.911 X 10^ 


4.677 X 10^(33.9) 
3.869 X 10^(0.028) 








70 < < 88 


52.8 


discrete 


fn^^ = 12 iV^^' 

V max 1 max 


= 24) 


0.758 


9.949 X 10^ 


7.541 X 10^ (54.7) 




full 




2 < < 17 


1.346 X 10^ 


1.527 X 10^ (11.1) 








n <E^.,< 70 


1.914 X 10^ 


4.683 X 10^(34.0) 
4.141 X 10^(0.30) 








70 < E^ < 190 


476.8 


discrete 


\ "max 1 max 


= 24) 


0.406 


1.850 X 10^ 


7.538 X 10^ (54.7) 
1.560 X 10^ (11.3) 




(l,r-r) 




2<E^< 17 


1.411 X 10^ 








17 < E^ < 70 


1.899 X 10^ 


4.653X W (33.7) 
4.140 X 10^(0.30) 








70 < E^ < 190 


475.6 



Table 3: The same quantities as in Table II, but for the operator (r^ — r]r)Y\^j, in ^'^^Pb. 
See the caption to Table II. 



Type of calculation 


energy in MeV 


B{D) in fm^ 


EWSR in fm^MeV (%) 


correct calculation 


0.0 

< 


0.0 


0.0 (0) 

5.63 X 10^ (40.8) 


continuum 
(1,t-t) 


0.107 

17 <E^< 70 
70 < < 150 


1.163 X 10^ 
8.294 X 10^ 
1.802 X 10^ 
0.005 X 10^ 


1.240 X 10^ (0.009) 
9.064 X 10^ (6.6) 

4.440 X 10<^ (32.2) 
0.042 X 10*^ (0.3) 


discrete = 12, TV^j; = 16) 
full 


1.401 
2 < < 17 
17 < < 70 
70 < < 88 


1.672 X 10^ 

9.369 X 10^ 
1.878 X 10-^ 
36.2 


2.342 X 10^ (0.011) 
1.006 X 10^ (7.3) 
4.580 X 10^ (33.2) 
2.653 X 10^ (0.019) 


discrete (nr. = 12, iV^j^ = 24) 
full 


0.758 
2 < < 17 
17 < < 70 
70 < < 190 


2.005 X 10^ 
9.492 X 10^ 
1.889 X 10^ 
368.1 


1.520 X 10^ (0.011) 
1.015 X 10^ (7.4) 
4.604 X 10^ (33.4) 
3.165 X 10^ (0.23) 


discrete (n^ = 12, iV^i^ = 24) 
(l,r.r) 


0.406 

2 < < 17 

17 < < 70 
70 < < 190 


3.607 X 10^ 
9.544 X 10^ 
1.886 X 10^ 
364.7 


1.464 X 10^ (0.011) 
1.023 X 10^ (7.4) 
4.597 X 10^ (33.3) 
3.130 X 10^ (0.23) 
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